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Signal Separation System and Method 

This invention relates to the field of signal separation, in particular to 
systems and methods for the separation of individual periodic source 
signals from a composite signal. 

5 Signals generated by many real-life sources can be characterised as 

periodic or almost periodic,. Typical examples include: sounds generated 
by the propellers of Ships, mechanical vibrations arising from rotating 
machinery and electrical signals as a result of a heart-beat. Such signals 
do not generally occur in isolation: multiple sources generating signals of 

10 the same nature (e.g. sound) are often present within the same 

environment and each will make its own contribution. Individual source 
signals may be similar 6r dissimilar in form but, regardless of this, 
discrimination between signals is likely to require more than straightforward 
monitoring of a detector response. The detector will respond to a 

15 composite signal which contains contributions from each source. In theory, 
the composite signal should be susceptible to mathematical analysis to 
retrieve its component parts. There are however many problems 
encountered in practice. 

Prior art analysis of a periodic composite signal is, in the first instance, 
20 dependent on estimation of the period of one or more of the signal 

components. The information that may potentially be available from an 
analysed signal has motivated the development of many period estimation 
techniques. The shape (amplitude structure) and period of a signal can 
provide useful information about its source. Applications have proved 
25 numerous: machine condition monitoring, target classification and 

communications, to name but a few. If each periodic sound in a composite 
signal generated by rotating machinery can be isolated and ascribed to a 
particular source, then the operational nature of the sources can be 
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individually monitored. Such monitoring can facilitate early diagnosis of a 
falling component, which in turn results in reduced downtime and repair 
costs. Detection and classification of a military target (e.g. a ship or 
submarine) is possible by analysis of the signals generated by its periodic 
5 emitters (e.g. propellers). Periodic modulation schemes are often used in 
telecommunications, and if multiple signals are transmitted on the same 
channel then analysis of the composite signal Is cleariy essential for 
decoding the individual messages. 

For all the above listed applications, it is important to be able to isolate 
10 individual periodic signals (which will be referred to herein as component 
signals) from a composite signal. Each individual signal potentially 
contains vital information relating to its source. 

One technique for separating such signals makes use of multiple sensors 
to detect the composite signal. Data collected within these multiple 
1 5 channels is then subject to an analysis known as Blind Signal Separation 
(BBS). This requirement for multiple sensors can be a major disadvantage: 
for separating n signals, at least n sensors are typically needed. Clearly 
the practicality of BSS is limited in many cases, with the limitation 
becoming more severe if a large number of source signals are present. 

20 The goal of BSS methods is the recovery (separation) of a number of 
statistically independent signals given only sensor obsen/ations that are 
mixtures of the source signals together with noise and unwanted 
interference. The term "blind" is used to indicate that no prior assumptions 
are made concerning either the nature of the sources (beyond their 

i 

25 statistical independence) or the manner in which their signals are 

combined. Description of BSS and some of its applications can be found in 
S. Haykin, "Unsupervised adaptive filtering, Volume 1 , Blind Source 
Separation", John Wiley and Sons, 2000 and A. Hyvarinen, J. Karhunen 
and E. Oja, "Independent Component Analysis", John Wiley and Sons, 



2001 . BSS makes use of the spatial characteristics (e.g. arrival direction) 
of each signal in order to separate component signals. 

Given the practical limitations of providing multiple, separated sensors, 
techniques have also been developed for separating a composite signal . 
detected at a single sensor position. Most prior art period estimation 
methods in fact exploit their potential for single-sensor analysis, making 
use of techniques like the Fourier transform. This is generally applied by 
the use of comb filters, which are tailored for application to periodic signals. 

Period estimation techniques aim, for each individual source signal present, 
first to provide an estimate of its period and then to extract an 
approximation to the shape of its base repeating unit. This base unit is the 
smallest repeating unit within the signal and it will be referred to herein as a 
"eyelet". If a signal is modulated by an amplitude envelope then the 
amplitude of each eyelet will vary as it repeats. Such modulations occur, 
for example, if a communication signal is sent on a earner wave or if a 
heart-beat is affected by exercise and each beat has a different strength. If 
multiple periodic signals are present in the composite data theii each 
individual signal is generally detected, estimated and subtracted in lum. 
until only noise remains. Most calculations are performed on a digitally 
sampled signal and hence the period may be measured in temis of the 
number of sample points that the eyelet spans. 

Further examples of modulating amplitudes can be found In: Y. Medan, 
E. Yair and D. Chazan, "Super Resolution Pitch Determination of Speech 
Signals", IEEE transactions on Signal Processing, Vol. 39, N0..1. January 
1991 in which periodic components of a speech signal are modelled with a 
varying scaling; M. Feder, "Parameter Estimation and Extraction of 
Helicopter Signals observed with a Wide-Band Interference", IEEE 
transactions on Signal Processing. Vol. 41 , N0.1 , January 1 993 in which 
helicopter acoustic signals are modelled with a varying sealing; and P. P. 



Kanjilal, S. Palit and G. Saha, "Fetal ECG Extraction from Single Channel 
Maternal ECG using Singular Value Decomposition", IEEE transactions on 
Biomedical Engineering, Vol. 44, N0;1, January 1997 in which heart-beat 
signals are modelled with a varying scaling. 

5 One technique used in the prior art to derive an estimate of the period of a 
signal is known as synchronous averaging. This technique Is based on the 
concept that period likelihood estimators can be used to measure the 
"strength" of period correlation for a range of period test values. It therefore 
requires an initial evaluation of a measure of likelihood for each of many 

10 test values of the period. The test value that results in the highest 

likelihood (or correlation) is taken to be a direct estimate of the true period. 
The signal is then partitioned in the time domain into eyelets of the 
estimated period and successive partitions averaged to obtain an estimate 
of their shape. Synchronous averaging is described by M. Feder in the 

15 reference cited above. 

The comb filter makes use of the principle that a periodic signal can be 
considered as being made up of basis periodic components at a 
fundamental frequency and multiples thereof. In the frequency domain, this 
method is equivalent to multiplying the Fourier transfomri of a signal with a 
20 filter that passes only the fundamental and harmonic frequencies of test 
periods. An example implementing this method is described by A. de 
Cheveigne and H. Kawahara in "Multiple period estimation and pitch 
perception model", Speech Communication, 27, 175 - 185, 1999. 

A third technique that may be used in the estimation process is singular 
25 value decomposition (SVD). Period SVD uses the same principle as 

synchronous averaging in that it first partitions a sampled signal in the time 
domain in accordance with a number of test periods. Instead of then 
averaging, an SVD is perfomied on the partitioned data. Partitioned 
segments are then combined using a weighted average and, at the correct 
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test period, will integrate coherently. Period SVD has the additional benefit 
that the modulation envelope (scaling of successive eyelets) of a periodic 
sigr^al is estimated along with the period. An example of its application is 
presented in the Kanjilal ef a/ reference, cited above. 

5 Many more examples of detection, estimation and separation of multiple 
periodic signals detected at a single sensor are found in the prior art. For 
example, see W.A. Sethares and T.W. Staley, "Periodicity transfonms". 
IEEE transactions on Signal Processing, 47(1), January 1999; V.P. 
Strokes, H. Lanshammar and A. Thorstensson, "Dominant pattern 

10 extraction from 3-D kinematic data", IEEE transactions on Biomedical 
Engineering, 46(1), January 1999. 

Every periodic signal can be considered as being made up of basis Fourier 
components at a fundamental frequency and at multiples of this 
fundamental (its harmonics). If periodic signals share common Fourier 
1 5 harmonics, this is tenned "overiap" In the frequency domain and signal 
leakage can occur. At the extreme, most significant overiap occurs if 
individual source signals within a composite signal have the same period. 

All the above-mentioned prior art decomposition techniques suffer from the 
fundamental problem that their accuracy is mari<edly decreased if 

20 component periodic signals overiap in the frequency domain. This arises 
because they are all based on using the signals' periods as the primary 
discriminant. For example, when using a comb filter, if there Is no 
frequency overiap then basis Fourier components at the fundamental 
frequency and harmonics will contain information relating to only one 

25 source signal. Combining the fundamental plus harmonics to reconstruct 
this source is therefore effective regardless of the presence or absence of 
other (different-frequency) signal components within the data. If overiap is 
present however then the process is not so clear cut. Very little distinction 
is available between basis Fourier components of different individual 




signals (originating from different sources), but which occur at the same 
position in the frequency domain. The principle of repeatedly subtracting 
individually detected and estimated signals Is no longer valid as each 
estimated signal will contain components (leakage) from other source 
5 signals. Such a scenario arises when the periods of component signals 
within the composite data are not coprime. The term "coprime" means that 
the highest common factoi^ between two positive integers (in this case the 
periods in sample point units) is one. That is, signal leakage will occur 
when the periods of one or more component source signals are, for 
10 example, the same or when one period is a multiple of the other. The 
problem is discussed by Z. Mouhan, C. Zhenming and R. Unbehauen. 
"Separation of periodic signals by using an algebraic method", IEEE 
ISCAS-91 , 5 pp 2427 - 2430, 1 991 . In this paper the authors highlight the 
inapplicability of their method to separating signals that are not coprime. 

15 It has to be borne in mind as well that as signals are digitally sampled for 
analysis, the signals need not be coprime only to the resolution of the 
sampling frequency for the prior art separation methods to break down. 
That is, signals of, for example, approximately equal periods, or 
approximately equal multiples, will also appear indistinguishable. 

20 This fundamental problem is not readily ignored as it is indeed the case 
that in many real-life situations, two or more component source signals will 
often occur with the same or similar periods. For example, if different parts 
of machinery rotate at the same speed then separate source signals will be 
generated at the rotation frequency. Similarly, if a signal reaches a sensor 

25 by multiple routes from its point of origin (multi-path), the overall received 
signal will be a combination of signals with the same period but having 
traversed different path lengths. Each path effectively gives rise to one 
perceived source and the combined signal is difficult to decompose into its 
individual components. Target classification is complicated in that the 

30 combined signal will not correspond to any known or actual signature. This 



problem is analysed in H. Amindvar and P.P. Moghaddam. "Estimation of^ 
propeller shaft rate and vessel classification in a multi-path environments". 
IEEE Proceedings of Sensor Array and Multi-channel Signal Processing 
Workshop, pp. 125 - 128. March 2000. In communications, optimal use of 
bandwidth is achieved if multiple signals are transmitted on the same . 
carrier frequency. Again the problem of separating individual messages 
with the same period is encountered. It would go some way to providing a 
complete process for signal analysis if, in the first instance, signals of equal 
or approximately equal periods received at a single sensor could be 
separated. 

In one attempt to avoid this type of problem a novel method of period 
estimation for communication signals was proposed by B. Santhanam and 
P. Maragos in "Harmonic analysis and restoration of separation methods 
for periodic signal mixtures: algebraic separation versus comb filtering", 
Signal Processing. 69 pp. 81 -91. 1998. Their approach exploits prior 
knowledge of the shapes of the periodic signals. Unfortunately, in most 
period estimation problems the shapes are unknown and so despite being 
found to enhance the detection and separation process, this method is of 
limited applicability. 

In summary therefore, when periodic signals significantly overlap in the 
frequency domain, the performance of existing period estimation 
techniques is degraded in two ways. Degradation is most severe when the 
signals have the same period. At a first stage of analysis, the ability .to 
even detect the presence of two signals is limited. If two signals having the 
same period are mixed, then only one may be detected using period 
discrimination alone. At a second stage, there is reduced capacity for 
signal separation. In general, the estimated separated signal will be a 
linear combination of the actual individual source signals. Current 
techniques to overcome this problem make use of an alternative 
discriminant: either arrival direction at. multiple sensors or the requirement 



of some foreknowledge of the shapes of the periodic component of the 
signals. Both of these techniques have significant disadvantages which 
prevent them from being universally applicable. 

There is a perceived need for a novel signal separation technique that is 
5 capable of separating individual source signals from composite data 

collected at a single sensor even in situations in which two or more source 
signals have very similar periods. 

Accordingly the present Invention provides a method of separating a 
plurality of source signals, each with similar or equal periodicity p, from a 
10 composite signal, the signal being in the form of a datastream with values 
indicative of sampled signal amplitude, the method comprising the steps of: 

(a) for a predetermined period p of the composite signal, forming a {/cx p) 
matrix X representing the datastream, wherein the sampled values are 
partitioned into consecutive segments of p sample points, and each 

1 5 row of X comprises the p values of one segment, data relating to 
consecutive segments being located in consecutive rows; 

(b) performing a singular value decomposition of the matrix Xto obtain two 
singular vector matrices U, 1/and a singular value matrix X\ and 

(c) perfonning an independent component analysis of one of the singular 
20 vector matrices U, VXo obtain an independent component matrix UR/, 

Ri^Vand an associated component matrix R^V, UkRi wherein one 
matrix UR2\ UkRi contains estimated separated signal envelopes af 
and the other /?AV^, /?/!/ contains estimated separated eyelets c/^ 

This method provides a technique by which a composite signal detected at 
25 a single sensor can be separated into its component individual source 
signals. This therefore avoids the requirement for multiple sensors and, 
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unlike prior art signal separation techniques, is specifically able to separate 
signals witli the same period from the composite. 

The predetermined period p is most preferably set at one of a number of 
test periods p'. steps (a) and (b) are repeated for each test period p', each 
5 derived matrix Xtest is assigned a likelihood indicator on the basis of values 
within its decomposed matrices Vand step (c) is only carried out on 
matrices derived from that matrix X which is assigned the maximum 
likelihood indicator. 

This method makes use of period SVD in an entirely novel signal 
1 0 separation problem. It provides the advantage of robustness to negative 
amplitudes in estimating a period within a composite signal, something that 
is not achieved by alternative prior art period estimation techniques. 

Alternatively, the predetermined period p may be estimated by 
synchronous averaging of the datastream. Although this method is not so 
1 5 robust as period SVD, it is more computationally efficient. 

in one embodiment of this invention, if the signal envelopes are more 
statistically independent than the eyelets then the independent analysis of 
step (c) is performed on the singular vector matrix Uto obtain the 
Independent component matrix t/A?/ containing estimated separated signal 
20 envelopes af and the matrix Rd^V containing estimated separated eyelets 

In an alternative embodiment, If the eyelets are more statistically 
independent than the signal envelopes then the independent analysis of 
step (c) is performed on the singular vector matrix VXo obtain the 
25 independent component matrix R/V containing estimated separated 
eyelets cf and the matrix UkRi containing estimated separated signal 
envelopes af. 



These two alternatives allow the method to be adapted to applications in 
which different aspects of the signals to be separated are statistically 
independent. 

The method may Include an additional step of estimating a number q of 
5 source signals of period p present within the composite signal and reducing 
the decomposed matrices £/, V, X in accordance with this number q of 
source signals prior to perfomning Step (c). This provides the advantage of 
reducing the computations necessary in carrying out the ICA stage of this 
invention. 

10 The number q of source signals may be estimated from a knowledge of the 
origins of these signals. Alternatively, it may be estimated from a number 
of elements within the singular value matrix X whose values exceed a 
predetermined threshold value. 

The composite signal is preferably detected at a single sensor. In an 
15 alternative embodiment the composite signal may be detected at at least 
two sensors. The matrix Xmay then, if mixing is instantaneous, be formed 
asa(Jkx p) matrix, where j is the number of sensors. Alternatively, one 
{kx p) matrix Xmay be fonned for each datastream. Each matrix X is then 
analysed independently and estimated source signals averaged. 

20 In an alternative aspect, this Invention provides a method of signal analysis 
comprising the steps of 

(a) separating a plurality of source signals into separated datastreams, 
each with a different period p; and 

(b) analysing each separated datastream by a method described 

25 hereinabove in order to determine the number of source signals with 

similar or equal periodicity within the datastream and, in the event that 
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this number is greater than 1 , to separate these source signals. 

In a third aspect this invention provides a method of machine condition 
monitoring comprising the following stages: 

(a) sensing a plurality of signals from a respective plurality of machinery 
5 sources, each signal with similar or equal periodicity p, at a sensor; 

(b) recording and digitising the sensor output to fomri a composite signal in 
the form of a datastream of values indicative of sampled signal 
amplitudes; 

(c) separating the plurality of source signals from the composite signal 
1 0 using a method as described above; and 

(d) analysing the separated source signals for indications as to the 
condition of respective machinery sources. 

In a further aspect the present invention provides a computer arranged to 
separate a plurality of source signals, each with similar or equal periodicity 
15 p, from a composite signal, the signal being in the form of a datastream 
with values indicative of sampled signal amplitude, the computer 
comprising: 

a datastream memory for storing the datastream; 

20 

a data partitioning module arranged to partition the datastream into 
representative (/c'x p) matrices X, wherein/for a number of values of a 
parameter p\ the sampled values are partitioned into consecutive 
segments of p' sample points, and each row of X comprises the p' values of 
25 one segment, data relating to consecutive segments being located in 
consecutive rows; 
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an SVD processing module arranged to perform a singular value 
decomposition of at least one of the matrices Xto obtain two singular 
vector matrices U, Vand a singular value matrix X; 

5 a period estimation module arranged to estimate the tme period p of the 
source signals from an average of data within rows of either the partitioned 
matrices Xor respective singular vector matrices K and 

an ICA processing module arranged to perform an Independent component 
1 0 • analysis of one of the singular vector matrices U, V generated by the SVD 
processing module (1 58) from the matrix X partitioned In accordance with 
the estimated period p and so to obtain an independent component matrix 
UR2\ fl/Vand an associated component matrix RzXV, UkRi wherein one 
matrix I//?/. UkRi contains estimated separated signal envelopes af and 
15 the other R^V, contains estimated separated eyelets cf. 

In a still further aspect the present Invention provides a computer-readable 
medium embodying instructions for execution by a processor, the 
instructions relating to separation of a plurality of source signals, each with 
similar or equal periodicity p. from a composite signal, the signal being in 
20 the form of a datastream with values indicative of sampled signal 
amplitude, the computer-readable medium comprising: 

(a) program code for. for a predetermined period p of the composite signal, 
forming a{kxp) matrix X representing the datastream, wherein the 
sampled values are partitioned into consecutive segments of p sample 

• 25 points, and each row of X comprises the p values of one segment, data 
relating to consecutive segments being located in consecutive rows; 

(b) program code for performing a singular value decomposition of the 
matrix Xto obtain two singular vector matrices U, Vand a singular 
value matrix X; and 
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(c) program code for performing an independent component analysis of 
one of the singular vector matrices U, Vto obtain an independent 
component matrix 1//?/. R/Vand ah associated component matrix 
R2XV, UkRi wherein one matrix i/f?/, UkRi contains estimated 
5 separated signal envelopes af and the other R^kV, fff'^V contains 

estimated separated eyelets cf. 

Embodiments of the invention will now be described by way of example 
only and with reference to the accompanying drawings. 

Figure 1 is a schematic illustration of an application of the method of this 
1 0 invention to machine condition monitoring. 

Figure 2a Is a schematic representation of a periodic signal. 

Figure 2b is a schematic representation of a periodic signal with varying 
modulations. 

Figure 3 is a representation of the fomiation of a mathematical model of a 
1 5 periodic signal with varying modulations. 

Figure 4 is a representation of the formation of a mathematical model of 
multiple periodic signals with varying modulations. 

Figure 5 is a flow chart Illustrating the steps involved in implementing a first 
embodiment of a method of this invention. 

20 Figure 6 is a representation of source eyelets used to generate data for an 
example implementation of a method of this invention. 

Figure 7a is a graphical plot of a single periodic signal. 

Figure 7b is a graphical plot of a composite signal containing three periodic 



signals, each of the same period and generated from the eyelets illustrated 
In Figure 6. 

Figure 8 is a plot of the eombination of the envelopes of the three periodic 
signals plotted in Figure 7. 

5 Figure 9 is. a plot of estimated source eyelets extracted from the data 

plotted In Figure 7 by a partial implementation of a method of this invention. 

Figure 10 illustrates plots of estimated amplitude envelopes extracted from 
the data plotted in Figure 7 by a partial implementation of a method of this 
invention. 

10 Figure 1 1 is a plot of estimated source eyelets extracted from the data 
plotted in Figure 7 by a full Implementation of a method of this invention. 

Figure 12 illustrates plots of estimated amplitude envelopes extracted from 
the data plotted in Figure 7 by a full implementation of a method of this 
invention. 

15 Figure 13 is a scatter plot illustrating graphically the processes carried out 
by a method of this invention in separating signals of the same or similar 
periods. 

Figure 14 is a schematic illustration of a computer system arranged to 
implement another aspect of this invention. 

20 Figure 15 is a representation of the formation of a mathematical model of 
multiple periodic signals with varying modulations detected at two sensors. 

Figure 16 is a series of four plots that illustrate the ability of a method of 
this invention to separate signals whose start times are misaligned. 
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Figure 1 shows a schematic illustration of the stages involved in a method 
of machine condition monitoring in accordance with an aspect of this 
invention. In the arrangement shown in this Figure, signals generated by 
three sources 1 0a. 1 0b, 1 0c of machinery rotating at the same speed are 

5 detected at a sensor 12. In one embodiment, the sensor 1 2 is a 

microphone that detects the sound of the rotating machinery 10a. b, c. In 
another, it is an accelerometer, which detects its vibration. In either case, 
the detected signal is converted by the sensor 12 (microphone or 
accelerometer) to an electrical signal. The precise nature of the signals 

10 themselves is not important. What is significant is that they originate from 
multiple sources and are periodic, or near periodic, with potentially more 
than one signal characterised by a given period. In the scenario shown in 
Figure 1 the structure of the recorded signal may be complex. Signals with 
the same period are present from different machinery components rotating 

1 5 at the same speed. 

As shown in Figure 1 , the process of machine condition monitoring can be 
broken down into four stages. The (composite) output signal from the 
sensor 12 is recorded and digitised at a data acquisition stage 14. This is 
followed by a signal separation stage 16 at which the composite signal is 

20 decomposed into its component parts. At an analysis stage 1 8 the 

individual signal components are analysed, for example for signs that the 
machinery 10a, 10b. 10c may have developed a fault. There are many 
existing analysis methods that can be applied to a single source signal to 
the benefit of machine condition monitoring. Examples are given by N. 

25 Arthur, J. Penman, A. McClean and A. Parsons in "Induction machine 
condition monitoring with higher order spectra I. Fundamentals and tixed 
frequency operation", IEEE 24* annual conference of the Industrial 
Electronics Society pp. 1889 - 1894. September 1998. Finally, if a fault is 
detected, corrective action is taken at a fault correcting stage 20. 
30 Depending on the fault, the source may be replaced or adjustments made 
to the machinery in order to reduce the operational stress on the source. 



It is to be noted that the first three stages above, namely data acquisition 
14, signal separation 16 and data analysis 18, are common to all problems 
involving multi-source (or composite) signal analysis. The novelty of this 
present invention lies in the signal separation stage 16 of analysis. The 
5 specific problem of machine condition monitoring is therefore used herein 
for illustrative purposes only and is not in any way to be considered limiting. 
This invention, in its broadest aspect, as will become apparent, has far 
wider application to any area of composite signal analysis. This is 
furthermore not limited to composite signals containing two or more . 

10 components at the same period; it is in this area however that 

improvements over the prior art will be most apparent. In other words, the 
novel technique of this present invention, which is capable of separating 
signal components of the same period can be applied more generally to a 
composite signal containing multiple components of the same and different 

1 5 periods. Prior art analysis techniques are first applied to extract the 
different-period components before application of a second analysis 
technique to each extracted component in order to separate same-period 
components hidden within. For clarity only therefore, it is assumed in this 
embodiment of the invention that different-period components are absent 

20 from the composite signal. 

Typical structures of periodic signals are illustrated schematically in 
Figure 2. In order to provide a mathematical basis for the signal processing 
carried out by this invention, it is necessary to set up a model for this 
structure. Figure 2a is an schematic representation of a perfect periodic 

25 signal 26 with period p. Four repeating units (eyelets) 28a, b, c, d, each 
extending through one period p are shown. The signal is digitised and so 
each eyelet 28 can be represented by p sample points, each having a 
particular amplitude Cn, n being an index parameter for the sampling. With 
reference to one eyelet 28a, this can be represented mathematically as a 

30 vector: 



c = 



C3 



Or, in shorthand, indicating a (p x 1 ) vector: (p x 1 ) c. 

A more general quasi periodic signal 30 is represented in Figure 2b. Again 
four eyelets 32a, b, c. d are shown, each spanning one period p of the 
signal 30. This signal 30 however is modulated by a time-varying 
envelope. As a result each eyelet 32a. b. e, d is scaled by a different 
amount: a first eyelet 32a by a factor a1 , a second 32b by a2, a third 32c 
by a3, a fourth 32d by a4, and so on. in effect, the modulation envelope 
can be considered sampled by each eyelet at its peak position. If the 
modulation envelope is sampled /f times within the segment that is 
analysed, a representation of the compJete signal can be generated by 
consideration of a block of k repeated eyelets. Modulation can be either 
smoothly-varying (as in amplitude modulation of a carrier wave) or 
discontinuous, by which it is meant that different and unrelated scaling 
factors are applied to successive eyelets (as, for example, in a heart-beat). 
In either case, the modulation envelope of a signal can provide vital 
information about its source. For example, in the situation described in 
relation to Figure 1 . if a piece of rotating machinery develops a fault, the 
fault may cause the machinery to come into contact with other machinery. 
Friction will increase and so more energy will be emitted by the source. 
This energy is in the form of an altered modulation envelope to the signal. 

Figure 3 is a representation of how an individual source signal is partitioned 
for processing by the method of this invention. In Figure 3a there is shown 

/c eyelets 32a. 32b 32k. each of period p. comprising a first source 

signal 34. The modulation envelope sampled at the peak of each eyelet 
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32a-k has amplitude aim, where the "1" index indicates that this is the first 
source signal under consideration and "rrf indexes the eyelet 32a-k within 
the block of eyelets 34. 

Figure 3b illustrates and indicates various fomis in which the block of k 
5 eyelets 34 can be represented as a matrix X. A graphical representation 36 
is shown towards the left hand side. The individual eyelets 32a-k within the 
block 34, each with their scaling factor a1 m, are partitioned and arranged 
vertically in a column 36. As there are p sample points within each eyelet, 
and /f scaling factors, Xis a (k x p) matrix. A mathematical representation 
10 38 has the scaling factors written as a column vector, each multiplied by c^, 
where denotes the transpose of the (p x 1 ) eyelet vector c. If a1_ is the 
(k x 1) vector representing the sampled points of the scaling envelope: 



all 
a12 

a^k 



then, mathematically, for a single source signal: 

15 X=a1.^ (1) 

As an aid to understanding, the (1 x p) vector c'' is represented graphically 
by the eyelet profile 32 at a third 40 representation of X. 

In Figure 4 the representation of Figure 3 is extended to two periodic 
signals, which have equal periods p. In Figure 4a, as before, the firet 

20 source signal 34 Is shown as /f eyelets 32a, 32b 32k, each of period 

p, and with amplitude aim, where "rrT indexes the eyelet 32a-k within the 
block of /f eyelets 34. A second source signal 42 Is also shown as ^eyelets 

44a, 44b , 44k, each of period p, but with a second profile 44. The 

modulation of this signal 42 is however different. Sampled at the peak of 
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each eyelet 44a-k the modulation envelope has amplitude a2m, where the 
"2" index indicates that this is the second source signal and "nf indexes the 
eyelet 44a-k within the block 42. 

Using the same partitioning model as before, the matrix X representing a 
composite signal consisting of both the first 34 and second 42 source 
signals is shown in Figure 4b. The scaling amplitudes a/y are now written 
into a (k X 2) matrix 46 and the two eyelet profiles (cj and one above 
the other, into a second matrix 48. Mathematically, the matrix X can be 
expressed as: 



X = AC = f^a,c] (2) 



where 3 Is a (k x 1) vector which contains the modulation amplitudes 
(sealings) of the fth eyelet and ^ is a (1 x p) vector which contains the 
eyelet profile of the /Ih periodic signal. Accordingly, A is a (k x 2) matrix 
containing the envelopes of the two periodic signals in its columns and C is 
15 a (2 X p) m,atrix containing the sampled eyelet profiles of each signal along 
its rows. 

When signals wrth the same period are present within a composite signal, 
their modulation envelopes will combine. The matrix A is therefore known 
as the mixing matrix as its coefficients operate on the source signals to 
20 form the received signal. 

Three distinct pieces of information are contained within the matrix X: the 
size of its rows (i.e. the period p of the source signals), the modulation 
envelopes (within A) and the eyelet profiles (within C). 

Figure 5 is a flow chart representation of the stages involved in separating 
25 signals in accordance with this invention. The separation process 1 6 will 
be referred to as the Blind Period Method (BPM). In applying BPM. it is 



) 

assumed that the periods of the source signals will be constant over the 
duration of the window used for detection and separation and also that 
synchronous sampling conditions will apply. That is, the periods will be 
integer-valued and the sampled points on each successive eyelet will be 
5 the same. It is to be noted that sampling rates may be varied to ensure 
that this latter condition applies. For example, a sample rate at which the 
period of a signal is measured at 89.25 samples can simply be increased 
by a factor of 4. This allows the period to be expressed as 357 samples at 
the new rate. 

10 The BPM essentially divides into three stages: a data partition stage 60, a 
Singular Value Decomposition (SVD) stage 62 and an Independent 
Component Analysis (ICA) stage 64. In a first embodiment, period 
estimation is carried out at the SVD stage 62. In further embodiments a 
different method Is used to carry out the period estimation prior to the SVD 

15 stage 62. In these embodiment the SVD stage 62 is only performed on 
data partitioned in accordance with this estimated period. 

In any embodiment of the analysis in accordance with the present 
invention, a datastream has been acquired in the form of signal amplitudes 
at a large number of sample points. It is assumed that multiple sources 
20 contribute to the signal but that they all have the same period p. 

With reference to Figure 5, in this first embodiment of the analysis method 
of this invention a first step 66 is, for each one of a set of test periods p\ to 
format the datastream In the form of a matrix Xtest- That is, for each value 
of the period p\ the first sampled amplitudes are written across the 
25 matrix Xtest in the first row. The next p' across the second row and so on 
until the matrix is complete at row k\ 

At step 68, SVD is carried out on the (k' x p*) matrix Xtest. SVD is a 
decorrelation and normalising method and is a well-known analytic 
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technique. The SVD of a general matrix X can be expressed as: 

X-UkV (3) 

where the (k x l<) matrix £/and the (k x p) matrix Vare orthononrial and X is 
a (k X k) diagonal matrix with real diagonal elements (the singular values) 
arranged in decreasing order. The columns of 1/ comprise the left 
orthonormal singular (temporal) vectors and the rows of ^comprise the 
right orthonormal singular vectors. Matrices t/and Vare unitary. 

It is important to understand what is achieved in carrying out an SVD on a 
matrix formed from the datastream of signal amplitudes. It is helpful in 
considering this to assume that the matrix X has been derived by correct 
partition of the datastream in accordance with the true period p of 
component signals present within. If partitioned in accordance with this 
period p it is assumed that the datastream fills k rows within matrix X The 
way in which the SVD of a set of matrices partitioned in accordance with a 
set of test periods leads to an indication of the true period will be explained 
later. 

In carrying out the SVD of the matrix X. the decomposition essentially 
obtains U, a matrix of amplitudes and V, a normalised waveform (each Is 
nomialised to unit power). In full, assuming k> p: 



20 X = 



^1 " 
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By multiplying out the first pair of matrices (t/ and X). it can be seen that the 
values within the /c columns are presen/ed, with the addition of a 
multiplicative factor o„, where n is the number of the column. Comparison 



of this (and equation (3)) with Figure 4b (and equation (2)) reveals that the 
columns of £7 above correspond to the columns of the matrix A in equation 
(2) and the rows of 1/ above correspond to the rows of the matrix C, That 
Is, infomiation relating to the eyelets {cj^ is contained in the rows of l/and 
information relating to the modulation envelopes (a^ is contained In the 
columns of U, The singular values within X indicate the strength of their 
interaction i.e. the strength of the estimated source signals. If the singular 
values are arranged in decreasing order, later columns in the Uk matrix will 
rapidly tend to zero - limiting the number of source signals found. In the 
event that q of the singular values are non-zero (g ^ /c), or can be 
considered as being above a threshold value, then U reduces to a (k x q) 
matrix, V denotes a (q x p) matrix and X denotes a (q x q) matrix. If q = 2, 
the comparison with Figure 4b is evident. 

Returning then to step 68 shown in Figure 5, it can be seen that, for each 
test period p\ performing SVD on the resulting matrix Xtest obtains sets of 
matrices £/, V^and X. Within the matrices V, each row comprises a 
weighted average of the relevant segment. At step 70 the rows within each 
matrix 1/are added together. At the correct test period p' = p, integration 
will be coherent, as the summation will inevitably include repeated eyelet 
profiles. If the test period is not equal to the true period, non-coherent 
integration results in a significantly lower average value. In this way a 
likelihood estimator is obtained for each value of the test period. The 
search commences at a low test period and increments systematically. 
Peaks will occur In a plot (referred to as a period likelihood estimate) of test 
period against segment average intensity when coherent integration 
occurs, that is, when the test period is equal to, or a'multiple of, the true 
period. An estimate of the true period p is therefore obtained at this step 
70. As a final step 72 of this period SVD stage 66, the decomposed 
matrices U, V and X that correspond to the matrix X partitioned in 
accordance with the true test period are extracted in preparation for the ICA 
stage 64 of the BPM. 
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It is not known in the prior art to apply SVD to the problem of separating 
signals with the same periods from a single-sensor datastream. In fact, it is 
only known that SVD may be applied to extract component signal periods 
from a composite signal. This feature is taught by the period SVD method 
of Kanjilal et ai, whose method is directed at signal separation when the 
periods may be coprime. The above-described partition 60 and SVD 62 
stages essentially follow this period SVD method but with the crucial 
difference that it is applied to data in which at least two component signals 
have the same period. This present invention is based on the realisation 
that if a single-channel datastream. representing values collected m the 
time domain, is partitioned in the manner described herein to fomi a matrix 
such as X, then the formerly unsolvable problem of separating signals of 
the same period collected in a single channel becomes solvable. 
Mathematically, it is in fact analogous to a problem that has been solved in 
the spatial domain, although far less intuitive to set up. That is, detecting a 
composite signal simultaneously at a number of detectors is more naturally 
written in matrix format, each row of the matrix corresponding to sampled 
composite signal values in a particular detector channel. For example, if 
two signals arrive at multiple sensors (at the same time) then the mixing 
model may be expressed as a combination of two amplitudes (scalings at 
each sensor) multiplied by the two signal components: 



Where the vector m; denotes the scalings of the /th signal at each sensor 
and denotes the shape of the /th signal. It is known that this problem 
(i.e. extracting /Wand Sfrom a matrix Xfomned from multi-channel 
datastreams) can be solved in the spatial domain by application of SVD 
and ICA (by BSS). it is not known that an analogous problem exists and 
can be solved in the time domain. The solution in a time domain is the 
subject of this invention. 
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As mentioned previously, it is l<nown in the prior art tliat a SVD of XUe 
matrix X will not generally reproduce the true signals in the spatial domain. 
Similarly it does not, in itself, reproduce the true signals in the time domain. 
This can be seen by reference to the example shown in Figures 6 to 10. 
5 Data within these examples has been simulated. It is representative of the 
situation illustrated in Figure 1 , but Is also generic to any application In 
which signals of the same period are mixed. 

Figure 6 is a plot 100 of three source eyelets Cj, Cg, cg, each of which 
belongs to a different periodic signal. The eyelets Cj, cg, Cg have the same 
10 period (40 sample points) and, although different, are highly correlated, as 
evidenced by their similar shapes. 

Figure 7a Is a plot of a section of a signal 102 based on one of the source 
eyelets Cj. Over a length of time the eyelet is replicated a number of 
times, with scaling in accordance with a signal modulation envelope. In this 
15 example, the modulation envelope corresponds to a real acoustic signal. 
The plot is generated by concatenating rows of the (k x 40) matrix 
X= [ai][c/], the modulation envelope being contained in the (k x 1) vector 
at and the eyelet profile in the (1 x 40) vector Cj^. 

A segment of a composite signal 104 generated by mixing three signals, 
20 each based on one of the source eyelets Cj, Cg, £3 Is plotted in Figure 7b. 
Within the time series shown, each eyelet cj, £g, £3 is replicated a number 
of times and source signals 102 are each modulated by a different 
envelope, which corresponds to a real acoustic signal such as speech or 
music. The composite signal 104 is generated by concatenating the rows 
25 of the (k x 40) matrix X: 



— I 



X = [a, 33] cl (5) 
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The (k X 1) vectors at. & contain the actual acoustic signals (speech 
and music) and the (40 x 1) vectors c^, Cg, £3 contain the eyelet profiles. 

Figure 8 is a pl6t of a composite acoustic signal 106. representing the 
mixing of the three envelopes . A single column of the composite matrix X 
5 represents a series of sample points taken on the composite envelope at 
the same posHlon on the eyelet. The composite acoustic signal 106 is 
accordingly generated by plotting values found in a single column of X as a 
time series. 

The SVD step 68 of this invention is then applied to the generated matrix X. 

1 0 Estimated eyelet profiles are contained within the rows of the matrix V and 
these results are shown in Figure 9. It can be seen that three eyelet 
profiles 108, 1 10, 112 are extracted from the decomposed matrix V. Not 
unexpectedly, as SVD applied in the spatial domain is known to be 
unsuccessful in separating signals of the same period, these eyelets 108. 

15 110.112 clearly do not correspond to the original source eyelets c^. cg, £3 
shown in Figure 6. That is, SVD applied in the time domain is also 
unsuccessful in separating signals of the same period. Modulation 
envelopes are contained within the columns of the matrix l/and these 
results are shown in Figure 10. Three modulation envelopes 114, 116. 118 

20 are estimated. From the similarities between two 1 1 6. 1 1 8 of the 

envelopes it is apparent that there is a large degree of signal leakage. 

The reason for the failure of SVD to separate the signals is that it 
constrains both the modulation envelopes and eyelet profiles to be 
uncorrelated. or dissimilar. This is inherent in the SVD methodology as it is 
25 a second order decorrelation method, which is intended to remove all 
similarities between signal pairs in a set of signals. Mathematically, this 
means that the decomposed vectors are made orthogonal. In many real 
signals, either the envelopes or eyelet profiles (as in this example) will be 
similar (correlated) and so a solution that makes them dissimilar will not 
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succeed in separating them. SVD also carries out a normalisation process, 
whicli forces normalised signals in the set to have the same power level. 
The singular values X reflect relative power levels and so allows their 
detection. It is important to note that, although the eyelet profiles are 
5 correlated, they do still need to be different for SVD to decompose the 
matrix X. 

SVD does however generally provide some indication of the number of 
source signals in the composite signal. At step 72, shown in Figure 5, it will 
be recalled that relative power levels of the component signals are 

10 contained in the singular values and these are generally arranged in 
decreasing order along the diagonal of matrix X. The number of source 
signals present, denoted by q, can be estimated by setting a threshold 
power level and determining the number of signals above it. In practical 
applications threshold estimation may be aided if the number of sources is 

15 known in advance. 

In summary therefore, period SVD, when applied to the problem of concern 
to this invention, first enables the correct period of component signals to be 
determined. Secondly, when the decomposed values of the correctly- 
partitioned matrix Xare retrieved, it enables the singular values to be used 
20 to determine whether more than one periodic signal with this period is 
present in the composite. 

Despite the inability of SVD to separate signals of the same period, it has 
been discovered that SVD does provide an initial basis for their separation. 
The SVD in fact produces linear combinations of the eyelets Cj, cg, cs, 
25 whichever combinations comprise orthogonal vectors. The goal of the third 
stage 64 of BPM is therefore to recover the actual eyelets and their 
envelopes from these combinations. 

Referring once again to Figure 5, the SVD stage 62, extracts the number q 
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Of sources and q linear combinations of eyelets from the composite signal. 
In some applications it is however difficult to set a suitable threshold level, 
and so it is assumed that ail singular values con-espond to sources. In 
such cases therefore either all decomposed signals are passed on to the 
next stage 64 or an estimated maximum number, for example, 30. To 
preserve the generality of the method it will be assumed in what follows 
that q,2^q<k signals are processed in the final stage 64. 

At the third stage 64. an Independent Component Analysis (ICA) is 
performed on the decomposed t/and V matrices. ICA is again an 
analytical technique that is known from BSS. BSS describes a process of 
extracting individual signals emitted by a number of statistically 
Independent sources from some linear combination received at an array of 
sensors. At the heart of ICA is its assumption that the combined sources 
are statistically Independent. To date. ICA has almost exclusively been 
applied to spatial problems: typically the separation of signals from multiple 
sources received at multiple sensors. A further assumption of JCA is that 
the signal combination mechanism Is linear, instantaneous and time- 
invariant. If these assumptions hold then the time-domain partitioning 
described above creates a mixing model that is analogous to the spatial- 
domain multiple sensor case. This third stage 64 of BPM is based on a 
novel, application of ICA to a time domain problem: a time sequence of 
data collected at a single sensor. In the prior art, ICA has not been applied 
to single-sensor period estimation problems. 

As ICA is based on an assumption of statistical independence, it is 
important to analyse a property of the composite signal that, as far as 
possible, exhibits statistical independence. A first step 76 of the ICA stage 
64 of BPM is therefore to determine whether it is the modulation envelope 
waveforms a, or eyelet profiles ci that are statistically Independent. In most 
applications it is envisaged that one skilled in the art of signal analysis will 
know this in advance. In the ease of machine condition monitoring it is 
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most likely that the envelopes are more independent: components of 
machinery rotating at the same speed are likely to have similar eyelet 
profiles, but what they connect or interact with (the source of the 
modulation) will depend on where a particular component is within the 
5 assembly. Similarly, if multiple speech signals are carried' in a single 
communications channel, the "acoustic" envelopes are unlikely to be 
correlated. On the other hand, in sonar applications, for example, in which 
a sonar pulse is corrupted by a multi-path signal, it is the eyelets that are 
more likely to be statistically independent. Since sonar can make use of 

1 0 narrow pulses then any time delay is likely to separate a pulse and its multi- 
path equivalent, which will therefore not overlap in the time domain. That 
is, the pulses will comprise eyelets with independent shapes in the time 
domain. Without independent knowledge of the shapes or delays, the 
signals of period p will appear to be one signal. In addition, the statistics of 

1 5 multi-path environments can change quickly and so only a few eyelet 
repeats would satisfy the requirement of having equal period for a valid 
separation. 

If it is the envelopes that are statistically independent, BPM then proceeds 
to step 78. This is to assume that the amplitudes, contained within matrix 
20 U, are related to the true independent envelopes by means of an unknown 
rotation Rz. Mathematically, X \s set at step 78 to: 

X= UR^R^y (6) . 

where R2 defines a unitary (q x q) rotation matrix and so RzRz = A the 
identity matrix. 

25 Once the matrix X has been fomiulated as described herein and then 
decomposed using SVD, at step 80 ICA is finally used to determine the 
hidden rotation matrix R2. As is well documented in relation to prior art 
applications of ICA, if signal combinations are linear. Instantaneous and 
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time-invariant, source signals cannot be separated using second order 
statistics. ICA is therefore based on higher order statistics (HOS), usually 
represented by cumulants of higher order than 2. 

ICA implemented in accordance with this embodiment of the invention is 
5 based on fourth order cumulants. Other implementations of ICA may also 
be used; it is not necessarily limited to fourth order statistics, but should 
make use of order three or greater. The procedure adopted to calculate R2 
from the matrix U is complex mathematically and well known in the prior 
art. Effectively R2 is computed as an iterative sequence of painwise 
1 0 rotations. Each rotation is designed to maximise the statistical 

independence of a given pair of modulation vectors a, within the matrix U. 
Statistical independence is measured by the corresponding painwise 
contrast function. The procedure adopted by this embodiment of the 
invention follows that described by J. G; McWhirter. I. J. Clarke and G. 
1 5 Spence in "Multilinear Algebra for Independent Component Analysis" 
SPIE's 44* Annual Meeting, the International Symposium on Optical 
Science, Engineering and Instmmentation. Denver. USA 1 8 - 23 July. 
1999. but other ICA algorithms can also be used. As mentioned previously 
ICA. Applied to spatial problems, is well known in the prior art. The specific 
20 details of applying ICA to the decomposed matrix X (as described herein) 
will therefore be obvious to one skilled in the art. The inventive aspect of 
the present invention lies in the realisation that, despite the prejudices of 
the prior art. SVD and ICA can be applied to a time domain problem in 
order to separate signals with the same period. 

25 Once the rotation rrjatrix R2 has been found, BPM computes, at step 82. 
final estimates of the eyelets c,* and modulation envelopes a,^ ^or q 
detected signals, the separated envelopes are contained In the columns of 
the (k X q) matrix t/ff/. The separated eyelet profiles are contained in the 
rows of the (q x p) matrix H.?. V. At step 84. the separated envelopes af 

30 and / or eyelets c;^ are passed to the analysis stage .1 8. 



~30~ 



In situations for which the modulation envelopes are not independent the 
ICA analysis based on U, as described above, breaks down. This is 
because the requirement for statistical independence is not met. As with 
SVD, ICA removes similarities within a dataset. If the modulations 
5 therefore are correlated, signals will not be correctly separated. In such 
cases, separation is still possible however, based on the independence of 
the eyelet profiles. Accordingly, if it is the eyelet profiles that are statistically 
independent, BPM proceeds to step 86, instead of step 78, after the 
selection step 76. The ICA has to be modified in order that analysis 
1 0 remains based on independent components - in this case the eyelet 
profiles. At step 86 matrix X is therefore set to: 

X= LfkRiR/V (7) 

where Ri defines a unitary (gx q) rotation matrix and so R-iR/^ I, the 
identity matrix. 

15 As it is the eyelet profiles as opposed to the modulation envelopes that are 
independent, Ri is calculated at step 88 from the matrix V. This calculation 
is carried out analogously to that described previously: Ri is computed as 
an iterative sequence of pairwise rotations, each rotation being designed to 
maximise the statistical independence of a given pair of eyelet vectors ci 

20 within the matrix V. 

Once the rotation matrix Ri is found, BPM computes, at step 90, final 
estimates of the eyelets cf and modulation envelopes a/^ For q detected 
signals, the separated envelopes are contained in the columns of the 
(k X q) matrix UkRi. The separated eyelet profiles are contained in the 
25 rows of the (q x p) matrix R/V. The alternative ICA processes then 
converge again and, at step 84, the separated envelopes af and / or 
eyelets cf are passed on to the data analysis stage 18- This stage 18 
employs known techniques that are used to extract information about the 



-31 ~ 



source of each signal component and to display results to a user. 

The full three stages of BPM were applied to the datastream used to 
generate the plots shown in Figures 6 to 10. Figure 1 1 displays a plot of 
three estimated eyelets cf, si, extracted from the composite signal 
5 1 04 by BPM (version that exploits the independence of the envelopes). 
Comparison with the source eyelets Ct, cg. of Figure 6. which were used 
to generate the composite signals in the first place, readily reveals the 
success of this method. 

IVIodulation envelopes are contained within the columns of matrix t//?/and 
1 0 the results extracted by BPM are shown in Figure 1 2. Three estimated 
modulation envelopes 122, 124, 126 are recovered. These acoustic 
signals are readily identifiable as having wavefonns characteristic of 
speech 122, 124 and music 126. In fad the con-elation coefficients 
between the orthononnal source modulations and their respective 
15 estimated modulations 122. 124, 126 are all greater than 0.99. The 
optimum correlation, representing perfect signal separation, is 1.0. 

The ICA stage 64 of BPM therefore proceeds by one of two alternative 
routes. First, if the eyelets are con-elated, signal separation (of eyelets and 
envelopes) is based on the independence of the envelopes. Secondly, if 
20 the envelopes are correlated, signal separation is based on the 

independence of the eyelet profiles. It is important to note that even 
although the eyelets or envelopes may be correlated, they must still be 
different for BPM to be effective. 

There is an additional advantage to be had In being able to use the 
25 independence of the eyelet profiles in the ICA analysis. Typically, to 

complete the analysis and make signals independent, at least 50 samples 
are required. If the independence of the envelopes is exploited, each, 
sample used in an envelope vector (element within a row of matrix U) is 



associated with one eyelet segment (period p samples long). Thus, for an 
envelope vector to span 50 samples, 50 replications of the eyelet will be 
required. By way of contrast, if the independence of the eyelets is 
exploited, signal separation may be achieved using far fewer eyelet 
5 repetitions. For q source signals, signal separation can be achieved using 
q eyelet segments (of period p, assuming that p > ~ 50) 

A helpful illustration of the mixing and then deeorrelation process is shown 
in Figure 13. Figures 13a-d are each referred to as a scatter plot in which a 
time-domain signal is sampled and an instantaneous amplitude (envelope) 

10 of one source signal is plotted against that of another, at each sample time. 
Thus if one source has amplitude 3 and the second has amplitude 1 at a 
particular sampling time, a point is placed on the plot at position (3,1). 
Over an extended interval, a number of points appear on the graph as the 
number of sample times increases. Figure 13a is a scatter plot 128 of two 

15 random source signals which are not mixed. As can be seen, points are 
distributed over the area of the plot 128. Figure 13b is a scatter plot 130 of 
a composite signal when the two source signals shown In Figure 13a are 
mixed. In this ease, data for one axis of the plot is taken from one column 
of the partitioned matrix X and data for the second axis is taken from 

20 another column. That is, the plot 130 shows a first combination of two 

envelopes, each weighted with one sample of their eyelets, plotted against 
a second combination of the envelopes. It can be seen that the plot 130 is 
far from an even distribution. Mathematically, the mixing in fact 
corresponds to a complicated combination of stretching, shearing and 

25 rotation. Figure 13c is a scatter plot 132 of two signals taken from the end 
of the SVD stage 62, 66 of processing. That is, the corresponding columns 
of Xto those used in Figure 3b are plotted after they have been made 
orthonormal: after X has been decomposed to a product UKV. It can be 
seen from this plot 132 that SVD removes the stretching and shearing 

30 caused by the mixing. Figure 13d shows a scatter plot 134 of the two 
columns after complete BPM processing. That is, of two columns of the 
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rotated matrix UR/. The effect of the hidden rotation n:iatrix found by the 
ICA stage 64 can be seen to be alignment of the straight edges of the 
Figure 3c scatter plot 132 with the co-ordinate axes. This representation 
applies regardless of whether the ICA stage is based on the independence 
of the eyelets or of the modulation envelopes. 

As an aid to summarising the signal separation process, components of a 
computer 1 50 arranged to carry out the separation and analysis are 
illustrated schematically in Figure 14. The computer 150 comprises an 
interface 152 to allow a user to select optional processing steps as 
required, a datastream buffer memory 154 for storing digitally sampled 
output from the sensor (the composite signal) and four processing modules 
156. 158, 159, 160 for perfonming the BPM separation. This embodiment 
of the invention acts on a signal with components of only one period p. It 
can however readily be adapted to perform a prior art period separation 
process before execution of the procedure described herein. The 
processing modules include a data partitioning module 156, with access to 
the datastream memory 154, which is arranged to partition data within the 
datastream memory 154 in accordance with a set of test periods ff; an 
SVD processing module 158 arranged to decompose the partitioned 
datastream into a number of linear combinations of source signals; a period 
estimation module 159 arranged to extract an estimate of the true period p 
of the periodic components of the composite signal and an ICA processing 
module 160 arranged to extract source envelope and eyelet profiles from 
the decomposed matrices £/. V, X obtained from the matrix X when it is 
partitioned in accordance with the true period p. In addition, the computer 
150 also includes two further memory units: a matrix memory unit 162 and 
a decomposed signal memory unit 164 for storing data output from the data 
partitioning module 156 and SVD processing module 158 respectively. In 
addition an analysis unit 1 66 is arranged to analyse the data for fault 
indication or other further processing function. 



In carrying out the BPM signal separation process as indicated In Figure 5, 
the sensor output in the form of sampled intensities is written to the 
datastream memory 154. The data partitioning module 156 accesses this 
data and, for each test period p', partitions the data to form a matrix Xtest- it 
5 then writes the data within each segment to a row of the matrix memory 
unit 162. The SVD processing module 158 accesses the data stored, as 
matrix Xtest, in the matrix memory unit 1 62 and perfomris a singular value 
decomposition of this matrix Xtest- Results, in the form of amplitude data, 
profile data and source strength data are written to the decomposed signal 

10 memory 164. For each matrix Xtest, the period estimation module 159 
calculates a weighted average of the profiles stored to obtain a likelihood 
Indicator of a particular test period p\ Once the entire set of test periods 
have been assigned a likelihood indicator, an estimate of the true period is 
made as that test period which has the maximum likelihood indicator. The 

1 5 decomposed data within the decomposed signal memory 1 64 that 

corresponds to that derived from the matrix X partitioned in accordance 
with the estimated period p is highlighted by the period estimation module 
159. The period estimation module 159 in addition accesses the source 
strength data in the decomposed signal memory 164 and estimates the 

20 number of signals that are present at this period p i.e. how many of these 
highlighted singular values exceed a preset threshold. The components of 
the decomposed matrices that correspond to the number of estimated 
sources are also highlighted in the decomposed signal memory 164. 

The ICA processing module 160 accesses the highlighted stored amplitude 
25 and profile elements for processing. Again either the amplitudes or profiles 
are processed, in response to a selection made by the user via the 
interface 152. After the ICA processing is complete, data relating to the 
separated eyelets and envelopes Is sent to the analysis unit 166 for 
analysis and / or display. 

30 In considering the BPM process described herein, it is apparent that If the 
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period were estimated in advance of the SVD stage 62 (Figure 5). then only 
a single SVD need be carried out on a con-ectly partitioned matrix. That is, 
we may be able to improve the computational efficiency of the analysis. As 
mentioned previously, period estimation techniques are well known in the 
5 prior art and any of a number of them can alternatively be used. Most of 
these period estimation techniques are based on the principle of 
synchronous averaging. A second embodiment (not shown) of this 
invention accordingly mal<es use of synchronous averaging. 

A synchronous averaging step is performed prior to the SVD stage 62. As 
1 0 in the previous embodiment, the datastream is divided into a number of 
segments, the size of each segment corresponding to a test period. If the 
test period is equal to the true period, coherent integration can be achieved 
simply by averaging the segments. If the test period is not equal to the true 
period, non-coherent integration results in a significantly lower average 
1 5 value. The test period which gives rise to the highest average value is 
therefore taken as an estimate of the true period p. This technique is 
described by Sethares et al. 

Thereafter, only the correctly-partitioned matrix X is passed to the SVD 
stage 62 of analysis. The ICA stage 64 is unchanged and proceeds on the 
20 basis of which parameter is statistically independent. 

The computer system 150 illustrated in Figure 14 is in fact capable of 
analysing the datastream in accordance with either embodiment of the 
invention. If the user selects period estimation by synchronous averaging, 
then the period estimation module 159 averages the data stored in the 
25 matrix memory 1 62 rather than in the decomposed signal memory 1 64. 
Note that this average is not weighted. Once the period is estimated, on 
the basis of the maximum likelihood indicator from the set of test periods p". 
the correctly-partitioned matrix is highlighted within the matrix memory 162. 
The SVD processing module 158 then perfomis an SVD of the highlighted 



matrix X and, as before, writes the output data to the decomposed signal 
memory 1 64. 

The disadvantage of synchronous averaging is that it cannot be used 
effectively if negative eyelets are anticipated. Cyclets with negative 
5 amplitudes (scalings) often arise in signals propagating through a medium. 
They are frequently encountered in sonar when multi-path reflections 
corrupt the original signal. In many other applications, the signal amplitude 
is measured about a subtracted mean and accordingly negative values 
result. In such cases, even at the true test period, positive and negative 

10 amplitude segments will average destructively. Coherent integration is 
limited and synchronous averaging may not generate distinctive peaks in 
the period likelihood estimate. By way of contrast, the period SVD method 
overcomes this limitation of synchronous averaging as it processes pairs of 
vectors; relative scalings of each eyelet pair are determined before they are 

15 combined. Accordingly the first-described embodiment that makes use of 
period SVD is preferred. 

As detailed above, the main benefit of BPM over prior art methods of signal 

separation is that it succeeds when applied to data received at only a single 

sensor. However this is not to say that multiple sensors cannot be used, 

20 indeed there is even a benefit to be had in reducing noise. In an 

implementation in which data is gathered at multiple sensors, BPM is 

applied in parallel to each datastream. The replicas of the estimated 

cyclets and / or envelopes obtained from each sensor recording may be 

averaged in order to suppress noise. For a eyelet repeated /c times across 

25 ./77 sensors, there are mk degrees of freedom to suppress noise, as 

I* 

opposed to only /c degrees in the single sensor case. There are however 
two distinct separation scenarios which must be considered before the 
averaging is done. 

In the first ease, referred to as instantaneous mixing, different sensor 
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outputs (channels) can be considered as consisting of replicas of the 
composite periodic signal but with different additive noise and different 
attenuations. The latter difference arises from the different propagation 
paths to spatially separate sensors. This scenario is represented In Figure 
15. In this Figure, two sources (not shown) generate one signal comprising 

eyelets ci and modulation amplitudes (a1 1 . a12, a1k) and another 

comprising eyelets Cg and modulation amplitudes (a21 , a22 a2k). 

Referring first to Figure 15a. there is shown the two source signals 170a, 
172a arriving at a first sensor 174 and the corresponding source signals 
170b, 172b, with different attenuation, arriving at a second sensor 176. 
The different attenuations can be represented by weighting the amplitudes 
of each signal 170a,b, 172a,b by a factor w/;, where /Indicates the source 
of origin and j indicates the sensor of reception. As Indicated In Figure 15b, 
In this case, segments from both 0' = 1 . 2, In this example) channels can 
simply be partitioned to form a {2k xp) (more generally (//f x p)) matrix X. 
BPM can then be carried out as before on the input matrix X. 

In the second scenario, the composite periodic signal arrives at different 
times at different sensors. Different sensor outputs can be considered as 
consisting of replicas of the periodic signal with different noise, attenuation 
and also time delays, the latter two both being caused by the different 
propagation paths to the sensor. Partitioning to use segments from all 
channels In the Input matrix X cannot be exploited in this ease. The time 
delays on different channels will mean that eyelet replicas from different 
channel recordings would be misaligned within the matrix X. • Attempting to 
separate such a misaligned matrix X has two consequences. First the 
number of source signals will appear to Increase. This arises as a delayed 
replica of a signal will appear to be distinct from the actual (undelayed) 
signal. Secondly, since replicas of the same eyelet are misaligned there 
will be non-coherent integration within the synchronous average step. In 
this scenario it is therefore necessary to perfonn two separate BPM 
analyses: one for each sensor channel. Time delays for each replica can 
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then be obtained using l<nown techniques and a con-ection applied. 
Corresponding source signals can then be averaged. 

In illustrating the analysis performed by the BPM, it has been assumed that 
the eyelets within each signal start at the same time. That is, the scaling 
5 factors a1 1 ,.a1 2, etc and a21 , a22, etc occur synchronously between 
respective signals. If this is not the case and the eyelets in a first signal 
begin at different times to those in a second signal, the BPM process is still 
able to separate the signals. As will be explained in more detail below, for 
each signal that is misaligned with the partitioning process, applying the 

10 BPM will result in an estimation of two separated signals. These two 
signals essentially correspond respectively to sum and difference 
components of the source sealing factors, with associated eyelets and 
singular values. In the majority of cases, most of the energy will be in the 
sum component; the difference component will be associated with a small 

1 5 singular value and can be neglected. Even If this Is not the ease, the BPM 
has only linearly decomposed the misaligned source signal and the two 
estimated (sum and difference) components may readily be associated and 
combined. 

To illustrate in further detail the BPM analysis of misaligned signals, 
20 consider first a situation in which two signals with the same period are 
mixed in a single time series. For this illustrative example, it is assumed 
that both replicated eyelets have a period of 100 samples but that the eyelet 
Cg has a different starting time from eyelet Cf . If we denote the samples 

within each eyelet vector c/'' as [c/(0), c,(1), c,<2) c/{99)], the amplitudes 

25 aij, as before, and the datastream itself as the (kp x 1 ) vector d, where the 
symbols kand p have the meanings assigned previously, then if cf is 
partitioned to fomn the matrix X, we have: 
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X = 



a1 1[o,(0). c,(1). c,(2) c,(99)] + a21Ic,(75).... c,(99)]a22[c,(0) 03(74)] 

a12[c,(0). c,(1), c,(2) c,(99)] + a22[c,(75).... c,(99)]a23[c,(0). ... c,i74)] 



(8) 



10 



where rt can be seen that each segment (of the two represented above) 
comprises the sampled amplitudes within the first eyelet c/. starting at 
t = 0 added to synchronously sampled amplitudes within the second eyelet 
c/. The second eyelet is displaced 25 sample points from the first. In this 
way. in the first segment, sampled points c^iTS) to 0^(99) are weighted by 
the amplitude a21 and added to the first 25 sannpled points of the first 
eyelet c/, whereas sampled points C2(0) to P2(74) are weighted by the 
amplitude a22 and added to the final 75 sampled points of the first eyelet 
c/. This is because the start of this eyelet c/ and the partitioning process 
are offset by 25 samples. 

From Equation (8) above it can be seen that each row of X contains a 
combination of the two eyelets and Xcan also be expressed as: 



15 



20 



X = 

321" 
a22 



an 

a12 



[Ci(0),c,(1),c,(2) c,(99)]+ 



[C2(75)....C2(99)A 0] + 



a22 
a23 



(9) 



[0, 0,C2(0),...C2(74)] 



where the second component has eyelet vector with 75 samples of zero 
amplitude and the third component has its first 25 samples of zero 
amplitude. Thus, the matrix X is effectively made up of three periodic 
components. The first periodic component consists of eyelet c, and the 
remaining periodic components consist of adjacent sections of the eyelet 
£2. This arrangement arises as a consequence of the starting point of £2 



not being aligned with tiie partitioning process and because modulations 
are varying, it can be obsen/ed that X is effectively a three rank matrix. 
Neither the amplitudes nor the eyelets of the latter two components (which 
belong to the same source signal) are independent: separation using BPM 
should, theoretically, prove difficult. 

If however the modulations of the eyelet £2 are slowly varying (i.e. if 
a21 « a22 « a23) then adjacent segments of £2 can be combined. The 
partitioned matrix Xean be expressed as: 



a11 
a21 



[c,(0),c,(1),c,(2). c,(99)] + 



a23 



[C2(75),...C2(99),C2(0)....C2(74)] 



(10) 



This formulation is similar to that shown in Equation (2) and the BPM can 
be used to separate the eyelets. Thus, effectively, with slowly varying 
modulations, the BPM will separate a non-aligned signal into one in which 
derived weightings are effectively averaged over the two eyelet amplitudes 
present in one segment. 



If, on the other hand, the modulation is highly non-stationary then the non- 
aligned signal is still separated, although not so straightforwardly, by the 
BPM. Such a highly non-stationary modulation can however be directly 
exploited by making use of the BPM version that is based on the 
independence of the modulation envelopes. 



Using the same notation as before, consider a single non-aligned signal 
with eyelets The partitioned matrix Xean be expressed as: 
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X = 



s2.\c^ (75), . . . Ca (99)]a22[C2 (0). . . . (74)]' 
a22[C2 (75).. . . (99)]a23[C2 (0).. . . (74)} 
a23[C2 (75), . . . (99)la24[C2 (0), . . . (74)] 



(11) 



In this example, X consists of two periodic components: the first example 
spanning samples 75 to 99 of £2 and the second spanning samples 0 to 74 
of £2. The modulations of two components in Equation (1 1) can be 
rewritten, by comparison with Figure 4 and Equation (2), as: 



A = [§2^3]= 



a21 a22 
a22 a23 
a23 a24 



(12) 



Where A denotes the mixing matrix. One of the modulation vectors 33 is 
delayed by one sample point with respect to the other modulation vector a^. 
If the modulations change abruptly then the two columns in Equation (12) 

10 will be sufficiently distinct to create two components in the analysed data, 
albeit that they will belong to the same replicated eyelet. Additionally, it is 
unlikely that these modulation vectors will be independent of each other. 
For this example, applying the BPM that is based on the statistical 
independence of the modulations, a pair of independent temporal 

15 modulation (basis) vectors would be generated that belong to the same 
signal. In order to recreate the time series for a single replicated eyelet, 
these independent modulation vectors will have to be combined. This is a 
straightforward task for a feingle signal. The BPM has simply linearly 
decomposed the data into a pair of independent temporal modulation 

20 vectors. The two signals can be expanded using the con-esponding 

singular values and the corresponding eyelets and simply added together. 

If a second signal is present then there is first a need to determine which of 
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the "independent" basis vectors are actually associated and so should be 
combined to recover the non-aligned signal. Despite this complication, the 
BPM has in fact achieved signal separation. It is only that two independent 
basis vectors are created for eyelets with a start that does not coincide with 
5 the partitioning process. 

One method of identifying and combining the associated base vectors may 
be appreciated by considering the following analysis. First, define the two 
separated (k x 1 ) modulation vectors of the same signal as 32 = a(k) and 
33 = a(k + 1), where the argument (k + 1) simply indicates that this vector is 
10 delayed by one sample with respect to a(k). Using a Taylor expansion up^ 
to first order for a about k, a(k + 1) can be expressed as: 

a(k+1)«a(k) + a'(k) (13) 

where a'(k) denotes the differential of a(k). The two vectors a(k) and 
a(k +1) can be more readily expressed as: 

[a(^+1) a(/f)]« [a{k) aik)+a\k)] 

[B{k)a{k) + a\k)Ma{k)a\k)]^^ ]j ^^"^^ 

where each modulation matrix is applied to modulations of the same eyelet. 
An alternative way to view the (k x 2) matrix consisting of a(k) and a(k + 1) 
Is the mixing of two components a(k) and a'(k) by a (2 x 2) matrix. In this 
ease, the BPIVI will estimate two modulation vectors and these are 

20 effectively the sum and difference components, albeit that they will be 
made independent of each other. Signaf association is though made 
possible by the fact that one modulation vector is effectively the 
(independent) differential of the other. This should enable their 
identification in the complete separated signals matrix. It is readily seen 

25 that if the modulations are slowly changing then most of the energy will be 
in the sum component and recombination of the sum and difference 
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components Is not required. That is. the situation reduces to the previous 
one in which an average amplitude was found. 

Figures 16a to d illustrate signal. separation by the BPM for a case in which 
the start times of the replicated eyelets (signals) differ. These plots (Figure 
16a - d) are simulated in a comparable way to the data shown In Figures 6 
to 12. Figure 16a is a plot 180 of three simulated source modulations 
182a. 182b, 182c. These are applied to respective replicated eyelets (not 
shown), which shall be referred to as cr, £g and £3. The offsets between 
the start of a signal eyelet cj, which is aligned with the partitioning process, 
and the starts of the other signals cg, ca are 1 1 and 27 samples 
respectively. Figure 16b is a plot 184 of the mixing of the three modulation 
envelopes. That is. It is generated by plotting values found in a column of 
matrix X. This Is in the same way as the plot 106 shown in Figure 8 is 
generated, but using the different modulations 182a, 182b, 182e and offset 
eyelets c,. cg, cg. After the SVD stage 62 is applied to the simulated matrix 
X. estimated modulation envelopes, contained within the columns of U, are 
shown in Figure 16c. Three modulation envelopes 188a, 188b, 188c are 
estimated. It is apparent that they 188a, 188b. 188c are not the same as 
the original source modulations 1 82a, 1 82b, 1 82c. By way of contrast, 
Figure 16d shows modulations 190a, 190b, 190c estimated using the BPM 
method. As stated above, the BPM estimates, for each of the offset 
signals, the largest (sum) component only. The modulations 190a. 190b, 
190c plotted in Figure 16d are evidently very similar to the original source 
modulations 182a, 182b, 182c shown In Figure 16a. In fact, the correlation 
coefficients between the orthonomnal source modulations and the 
estimated "sum" modulations were measured at 1 , 0.93 and 0.95. It can 
thus be seen that there is very little enror introduced by estimating only the 
"sum" component of a non-aligned signal in the BPM analysis. 
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Claims 

1 . A method of separating a plurality of source signals, each with similar 
or equal periodicity p, from a composite signal (104), the signal being in 
the form of a datastream with values indicative of sampled signal 
amplitude, the method comprising the steps of: 

(a) for a predetenmined period p of the composite signal, forming a 

(/f X p) matrix X representing the datastream, wherein the sampled 
values are partitioned into consecutive segments of p sample 
points, and each row of X comprises the p values of one segment, 
data relating to consecutive segments being located in consecutive 
rows; 

(b) performing (68) a singular value decomposition of the matrix Xto 
obtain two singular vector matrices U, l^and a singular value matrix 
"k; and 

(c) perfomning (80, 88) an independent component analysis of one of 
the singular vector matrices U, V\o obtain an independent 
component matrix UR/, R/Vand an associated component matrix 
Rz^V, UkRi wherein one matrix UR/, UkRi contains estimated 
separated signal envelopes and the other R2XV, contains 
estimated separated eyelets cf. 

2. A method according to claim 1 wherein the predetermined period p is 
set at one of a number of test periods p', steps (a) and (b) are repeated 
for each test period p' , each derived matrix Xtest is assigned a 
likelihood indicator on the basis of values within its decomposed 
matrices V, X and step (c) Is only carried out on matrices derived from 
that matrix X which is assigned the maximum likelihood Indicator. 
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3. A method according to claim 1 wherein the predetennined period p is 
estimated (68) by synchronous averaging of the datastream. 

4. A method according to claim 1 , 2 or 3 wherein the signal envelopes are 
more statistically independent that the eyelets and the independent 
analysis of step (c) is performed (80) on the singular vector matrix (/to 
obtain the independent component matrix y/?/ containing estimated 
separated signal envelopes ^ and the matrix B^XV containing 
estimated separated eyelets,©*. 

5. A method according to claim 1 , 2 or 3 wherein the eyelets are more 
statistically independent that the signal envelopes and the Independent 
analysis of step (c) is perfomied (88) on the singular vector matrix Vto 
obtain the independent component matrix f?/V containing estimated 
separated eyelets si and the matrix £/^/?t containing estimated 
separated signal envelopes a/". 

6. A method according to any preceding claim wherein the method 
includes an additional step of estimating a number q of source signals 
of period p present within the composite signal (104) and reducing the 
decomposed matrices t/. V, X in accordance with this number q of 
source signals prior to performing Step (e). 

7. A method according to claim 6 wherein the number q of source signals 
Is estimated from a knowledge of the origins of these signals. 

8. A method according to claim 6 wherein the number q of source signals 
is estimated from a number of elements within the singular value matrix 
X whose values exceed a predetermined threshold value. 



9. A method according to any preceding claims wherein the composite 
signal (104) Is detected at a single sensor (12). 



1 0. A method according to any one of claims 1 to 8 wherein the composite 
signal (104) is detected at least two sensors (174, 176) and the matrix 
Xls formed as a (//ex p) matrix, where j is the number of sensors. 

If 

1 1 . A method according to any one of claims 1 to 8 wherein the composite 
signal (104) is detected at least two sensors (174, 176) and the matrix 
Xis analysed for the datastream collected at each sensor (174, 176). 

12. A method of signal analysis comprising the steps of 

(a) separating a plurality of source signals into separated datastreams, 
each with a different period p; and 

(b) analysing each separated datastream by the method of any 
preceding claim in order to determine the number of source signals 
with similar or equal periodicity within that datastream and, In the 
event that this number is greater than 1 , to separate these source 
signals. 

13. A method of machine condition monitoring comprising the following 
stages: 

(a) sensing a plurality of signals from a respective plurality of 
machinery sources, each signal with similar or equal periodicity p, 
at a sensor (12); 

(b) recording and digitising the sensor (12) output to form a composite 
signal in the form of a datastream of values indicative of sampled 
signal amplitudes; 



(c) separating the plurality of source signals from the composite signal 
using a method in accordance with any one of claims 1 to 1 1; and 
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(d) analysing the separated source signals for indications as to the 
condition of respective machinery sources. 

14. A computer (150) arranged to separate a plurality of source signals, 
each with similar or equal periodicity p, from a composite signal (104), 
the signal being in the form of a datastream with values indicative of 
sampled signal amplitude, the computer (150) comprising: 

a datastream memory (154) for storing the datastream; 

a data partitioning module (156) arranged to partition the datastream 
into representative (/c'x p) matrices X. wherein, for a number of values 
of a parameter p', the sampled values are partitioned into consecutive 
segments of p' sample points, and each row of X comprises the p' 
values of one segment, data relating to consecutive segments being 
located in consecutive rows; 

an SVD processing module (158) arranged to perfomn a singular value 
decomposition of at least one of the matrices Xto obtain two singular 
vector matrices U, Vand a singular value matrix X; 

a period estimation module (159) arranged to estimate the true period p 
of the source signals from an average of data within rows of either the 
partitioned matrices Xor respective singular vector matrices V; and 

an ICA processing module (160) arranged to perfomi an independent 
. . component analysis of one of the singular vector matrices U, V 
generated by the SVD processing module (158) from the matrix X 
partitioned in accordance with the estimated period p and so to obtain 
an independent component matrix t//?/, R/Vand an associated 
component matrix RzXV, UkRi wherein one matrix l/fl/. UkRi contains 
estimated separated signal envelopes af and the other RzXV, R/V 



contains estimated separated eyelets c/^ 

15. A computer-readable medium embodying Instructions for execution by 
a processor, the instructions relating to separation of a plurality of 
source signals, each with similar or equal periodicity p, from a 
composite signal (104), the signal being in the form of a datastream 
with values indicative of sampled signal amplitude, the computer- 
readable medium comprising: 

(a) program code for, for a predetermined period p of the composite 
signal, forming a (/cx p) matrix X representing the datastream, 
wherein the sampled values are partitioned into consecutive 
segments of p sample points, and each row of X comprises the p 
values of one segment, data relating to consecutive segments 
being located in consecutive rows; 

(b) program code for performing a singular value decomposition of the 
matrix Xto obtain two singular vector matrices t/, l^and a singular 
value matrix X; and 

(c) program code for performing an independent component analysis 
of one of the singular vector matrices t/, Vto obtain an 
independent component matrix (/ffe^, R/Vand an associated 
component matrix R2XV, UKR^ wherein one matrix t//?/, UKRi 
contains estimated separated signal envelopes af and the other 
R2XV, /?7^1/ contains estimated separated eyelets cf. 
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